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INTRODUCTION 

Aquatic  invertebrates  are  aptly  applied  to  bioassessment  since  they  are  knoMvn 
to  be  important  indicators  of  stream  ecosystem  health  (Hynes  1970).  Long  lives, 
complex  life  cycles  eind  limited  mobility  mean  that  there  is  ample  time  for  the  benthic 
conxmvmity  to  respond  to  cximulative  effects  of  enviroimiental  perturbations.  This  report 
summarizes  data  collected  in  2003  from  4  sites  on  Bear  Creek  in  the  Gallatin 
watershed.   Sites  were  located  in  Gallatin  County,  Montana  and  lie  within  the  Middle 
Rockies  ecoregion  (Woods  et  al.  1999). 

A  multimetric  approach  to  bioassessment  such  as  the  one  appHed  in  this  study 
uses  attributes  of  the  assemblage  in  an  integrated  way  to  measxire  biotic  health.  A 
stream  with  good  biotic  health  is  "...a  balanced,  integrated,  adaptive  system  having  the 
full  range  of  elements  and  processes  that  are  expected  in  the  region's  natural 
environment..."  (Karr  and  Chu  1999).  The  approach  designed  by  PlaQdn  et  al.  (1989) 
and  adapted  for  use  in  the  State  of  Montana  has  been  defined  as  ". . .  an  array  of 
measures  or  metrics  that  individually  provide  information  on  diverse  biological 
attributes,  and  when  integrated,  provide  an  overall  indication  of  biological  condition." 
(Barbour  et  al.  1995).  Community  attributes  that  can  contribute  meaningfully  to 
interpretation  of  benthic  data  include  assemblage  structure,  sensitivity  of  cormnunity 
members  to  stress  or  pollution,  and  fianctional  traits.  Each  metric  component 
contributes  an  independent  measure  of  the  biotic  integrity  of  a  stream  site;  combining 
the  components  into  a  total  score  reduces  variance  and  increases  precision  of  the 
assessment  (Fore  et  al.  1996).  Effectiveness  of  the  integrated  metrics  depends  on  the 
applicability  of  the  vuiderlying  model,  which  rests  on  a  foundation  of  three  essential 
elements  (BoUman  1998a).  The  first  of  these  is  an  appropriate  stratification  or 
classification  of  stream  sites,  typically  by  ecoregion.  Second,  metrics  must  be  selected 
based  upon  their  ability  to  accxirately  express  biological  condition.  Third,  an  adequate 
assessment  of  habitat  conditions  at  each  site  to  be  studied  enhances  the  interpretation 
of  metric  outcomes.  '^' 

ImpUcit  in  the  multimetric  method  and  its  associated  habitat  assessment  is  an 
assumption  of  correlative  relationships  between  habitat  measvires  and  the  biotic 
metrics,  in  the  absence  of  water  quality  impairment.  These  relationships  may  vary 
regionally,  requiring  an  examination  of  habitat  assessment  elements  and  biotic  metrics 
and  a  test  of  the  presumed  relationship  between  them.  BoUman  ( 1998a)  studied  the 
Eissemblages  of  the  Montana  Valley  and  Foothill  Prairies  (MVFP)  ecoregion  and  has 
recommended  a  battery  of  metrics  applicable  to  the  montane  ecoregions  of  western 
Montana.  This  metric  battery  has  been  shown  to  be  sensitive  to  impairment,  related  to 
m,easures  of  habitat  int^rity,  and  consistent  over  repHcated  samples.  However,  scoring 
criteria  developed  for  the  MVFP  ecoregion  metric  battery  miay  not  be  appropriately 
sensitive  for  streams  of  the  Middle  Rockies  ecoregion.  Additional  work  may  be  needed  to 
calibrate  scoring  to  montane  regions,  thus  impairment  classifications  must  be 
interpreted  with  care.   For  comparison,  the  bioassessment  index  found  in  the  Montana 
Department  of  Environmental  Quality  (Montana  DEQ)  Standard  Ojjerating  Procedvires 
(Bukantis  1998)  is  also  scored  and  presented.  However,  this  index  has  not  been  tested 
for  sensitivity  to  impairment,  relationship  to  measures  of  habitat  integrity,  or 
consistency  over  repUcated  samples. 

In  this  report,  3  assessment  methods  are  used:  first,  taxonomic  data  is 
translated  in  to  2  bioassessment  indices,  and  m.etric  scores  in  each  index  are  siunmed 
to  derive  impairment  classifications  and  use  support  designations.  Second,  a  narrative 
interpretation,  based  on  the  author's  professional  judgment  is  given.  Metric  , 

performance  and  taxonomic  data  are  both  applied  to  this  analysis.  Third,  the  model  of 
Barbour  and  Stribling  (1991)  is  applied  to  bioassessment  and  habitat  assessment 
scores.  While  the  interdependence  of  these  methods  is  obvious,  since  the  same  data  are 
used  for  all,  some  degree  of  independence  is  maintained  throughout  the  analysis. 
Narrative  interpretations  are  given  without  regard  to  the  hioassessment  index  result, 
and  without  reference  to  habitat  assessment.  Since  indices  are  sunmiations,  they  can 


often  mask  evidence  of  impairment;  the  narratives  attempt  to  expose  the  potential 
shortcomings  of  the  indices.  Similarly,  graphing  the  eussociation  between  habitat 
assessment  scores  and  bioassessment  scores  using  the  model  of  Barbour  and  Stribling 
can  provide  clues  to  offer  support  or  possible  refutation  of  the  conclusions  of  the 
narrative  anafysis. 

BSBTHODS 

Samples  were  collected  in  July  2003  by  personnel  of  the  Gallatin  Local  Water 
Quality  District.  Sample  designations  and  site  locations  are  indicated  in  Table  1  and  on 
Figure  1 .  The  site  selection  and  kicknet  sampling  method  employed  were  those 
recommended  in  the  Montana  DEQ  Standard  Operating  Procedures  for  Aquatic 
Macroinvertebrate  Sampling  (Bukantis  1998).  Aquatic  invertebrate  samples  were 
delivered  to  Rhithron  Associates,  Inc.,  Missoula,  Montana,  for  laboratory  and  data 
.    analyses.  In  the  laboratory,  the  Montana  DEQ-recommended  sorting  method  was  used 
to  obtain  subsamples  of  at  least  300  org£uiismis  &om  each  sample,  when  possible. 
Organisms  were  identified  to  the  lowest  possible  taxonomic  levels  consistent  with 
Montana  DEQ  protocols. 

pr  Table  1.  Sample  designations  and  locations.  Sites  are  listed  in  upstream- to- 
downstream  order.  Bear  Creek.  2003. 


Site  ID 


Station  n> 


Activity  ID         Location  Descrtption 


Collection 
Date 


LaUtude 
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04 
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BEAR  CREEK  50  FT 
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BEASC 
05 


BEAR  CREEK  ACROSS 
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110*-55'-2  3.624" 
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BEAR  CREEK  0.25  MI 
M05BEARC03       03-U349-M        ABOVE  BOZEMAN  TAIL        07/24/03       45'-39'-3.025" 

ROAD 


110' -57" -9. 336" 


Two  bioassessment  indices  were  employed  in  evaluating  the  data.  The  first  index 
employed  was  developed  for  streams  of  western  Montana  ecoregions  (Bollman  1998a). 
Each  metric  used  in  this  index  was  tested  for  its  response  or  sensitivity  to  varying 
degrees  of  hvunan  influence.  Correlations  have  been  demonstrated  between  the  metrics 
and  various  symptoms  of  humam-caused  impairment  as  expressed  in  water  quality 
parameters  or  instream,  streambank,  and  stream  reach  morphologic  features.  Metrics 
were  screened  to  minimize  variability  over  natural  environmental  gradients,  such  as  site 
elevation  or  sampling  season,  which  might  confound  interpretation  of  results  (Bollman 
1998^.  The  multimetric  index  incorporates  multiple  attributes  of  the  sampled 
assemblage  into  an  integrated  score  that  accurately  describes  the  benthic  community  of 
each  site  in  terms  of  its  biologic  integrity.  In  addition  to  the  metrics  included  in  the 
index,  other  metrics  shown  to  be  applicable  to  biomonitoring  in  other  regions  (Kleindl 
1995,  Patterson  1996,  Rossano  1995)  were  used  for  descriptive  interpretation  of  results. 
These  metrics  include  the  number  of  "dinger"  taxa,  long-lived  tcixa  richness,  the  percent 
of  predatory  organisms,  and  others.  They  are  not  included  in  the  integrated 
bioassessment  score,  however,  since  their  perfonnance  in  western  Montana  ecoregions 
is  unknown.  However,  the  relationship  of  these  metrics  to  habitat  conditions  is  intuitive 
and  reasonable. 
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Figure  1.  Approximate  locations  of  sampling  sites.  Bear  Creek,  2003. 


The  six  metrics  included  in  the  bioassessment  index  used  for  MVFP  sites  in  this 
stucfy  were  selected  because,  both  individually  and  as  an  integrated  metric  battery,  they 
are  robust  at  distinguishing  impaired  sites  from  relatively  unimpaired  sites  (Bollman 
1998a).  In  addition,  they  are  relevant  to  the  kinds  of  impacts  that  are  present  in  the 
Upper  Kootenai  drainage.  They  have  been  demonstrated  to  be  more  variable  with 
anthropogenic  disturbance  than  with  natural  environmental  gradients  (Bollman  1998s^. 
Each  of  the  six  metrics  developed  and  tested  for  western  Montana  ecor^ons  is 
described  below. 

1 

1.  ESphemeroptera  (mayfly)  taxa  richness.  The  number  of  mayfly  taxa 
declines  as  water  quality  diminishes.  Impairments  to  water  quality  which  have  been 
demonstrated  to  adversely  affect  the  ability  of  mayflies  to  flourish  include  elevated 
water  temp>eratures,  heavy  metal  contamination,  increased  turbidity,  low  or  high 
pH,  elevated  specific  conductance  and  toxic  chemicals.  Few  mayfly  species  are  able 
to  tolerate  certain  disturbances  to  instream  habitat,  such  as  excessive  sediment 

■  deposition. 

2.  Plecoptera  (stonefly)  tajca  richness.  Stoneflies  are  particvdarly 
susceptible  to  imipairments  that  affect  a  stream  on  a  reach-level  scale,  such  as  loss 

'   :     of  riparian  canopy,  streambank  instability,  channelization,  eind  alteration  of 

morphological  features  such  as  pool  frequency  and  function,  rifQe  development  and 
sinuosity.  Just  as  all  benthic  organisms,  they  are  also  susceptible  to  smaller  scale 
habitat  loss,  such  as  by  sediment  deposition,  loss  of  interstitial  spaces  between 
substrate  particles,  or  unstEible  substrate. 

3.  Trichoptera  (caddisfly)  taxa  richness.  Caddisfly  taxa  richness  h£is  been 
shown  to  decline  vrtien  sediment  deposition  affects  habitat.  In  addition,  the  presence 
of  certain  case-building  caddisflies  can  indicate  good  retention  of  woody  debris  and 
lack  of  scouring  flow  conditions. 

4.  Number  of  sensitive  taxa.  Sensitive  taxa  are  generally  the  first  to 
disappear  as  anthropogenic  distvirbances  increase.  The  list  of  sensitive  taxa  used 
here  includes  organisms  sensitive  to  a  wide  range  of  disturbances,  including  warmer 
water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment 
deposition,  substrate  instability  and  others.  Unimpaired  streams  of  western 
Montana  typically  support  at  least  four  sensitive  taxa  (Bollman  1998a). 

5.  Percent  filter  feeders.   Filter-feeding  organisms  are  a  diverse  group; 
they  capture  small  particles  of  organic  matter,  or  organically  enriched  sediment 
xaaterial,  from  the  water  column  by  means  of  a  variety  of  adaptations,  such  as 
silken  nets  or  hairy  appendages.  In  forested  montane  streams,  filterers  are  expected 
to  occur  in  insignificant  numbers.  Their  abundance  increases  when  canopy  cover  is 
lost  and  when  water  temperatures  increase  and  the  accompanying  growth  of 
filamentous  algae  occurs.  Some  filtering  organisms,  specifically  the  Arctopsychid 
caddisflies  (Arctopsyche  spp.  and  Parapsyche  spp.)  build  silken  nets  with  large  mesh 
sizes  that  capture  small  organisms  such  as  chironomids  and  early-instar  mayflies. 
Here  they  are  considered  predators,  and,  in  this  study,  their  abundance  does  not 
contribute  to  the  percent  filter  feeders  metric. 

6.  Percent  tolerant  taxa.  Tolerant  taxa  are  ubiquitous  in  stream  sites,  but 
\rtien  disturbance  increases,  their  abundance  increases  proportionately.  The  Ust  of 
taxa  used  here  includes  organisms  tolerant  of  a  wide  range  of  disturbances, 
including  warmer  water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution, 
sediment  deposition,  substrate  instability  and  others. 

The  second  index  used  here  is  the  metric  battery  recommended  in  the  Montana 
DEQ  standard  operating  procedures  (Bukantis  1998).  These  metrics  should  be 
considered  provisional,  since  correlative  relationships  between  them  and  meaningful 
measures  of  habitat  condition  and  water  quality  have  not  been  evaluated.  Thus, 


conclusions  concerning  bioassessment  based  upon  these  metrics  must  be  regarded  as 
tentative. 

Scoring  criteria  for  each  of  the  metrics  are  presented  in  Tables  2a  and  2b.  Table  2a 
gives  metrics  and  criteria  for  the  revised  bioassessment  index  (Bollman  1998a)  and 
Table  2b  gives  metrics  and  criteria  for  the  Montana  DEQ  method  (Bukantis  1998). 
Metrics  differ  in  their  possible  value  ranges  as  well  as  in  the  direction  the  values  move 
as  biological  conditions  change.  For  example,  Ephemeroptera  richness  values  may 
range  from  zero  to  ten  taxa  or  higher.  Larger  values  generally  indicate  favorable  biotic 
conditions.  On  the  other  hand,  the  percent  filterers  metric  may  range  from  0%  to  100%; 
in  this  case,  larger  values  are  negative  indicators  of  biotic  health.  To  facilitate  scoring, 
therefore,  metric  values  were  transformed  into  a  single  scale.  The  range  of  each  metric 
has  been  divided  into  four  parts  and  assigned  a  point  score  between  zero  and  three.  A 
score  of  three  indicates  a  metric  value  similar  to  one  characteristic  of  a  non-impaired 
condition.  A  score  of  zero  indicates  strong  deviation  from  non- impaired  condition  and 
suggests  severe  degradation  of  biotic  health.  Scores  for  each  metric  were  summed  to 
give  an  overall  score,  the  total  bioassessment  score,  for  each  site  in  each  sampling 
event.  These  scores  were  expressed  as  the  percent  of  the  maximum  possible  score.  The 
total  bioassessment  score  for  each  site  was  expressed  in  terms  of  use- support.  Criteria 
for  use-support  designations  were  developed  by  Montana  DEQ  and  are  presented  in 
Table  3  a.  Scores  were  also  translated  into  impairment  classifications  according  to 
criteria  outlined  in  Table  3b.  * 

Table  2a.  Metrics  and  scoring  criteria  for  bioassessment  of  montane  streams  of  Western 

Montana  (Bollman  1998a). 

Score 

Metric  3  2  1  0 


Ephemeroptera  taxa  richness 

>5 

5  -4 

3-Q 

\ 

<2 

Plecoptera  taxa  richness 

>3 

3-2 

1 

0 

Trichoptera  taxa  richness 

>4 

4-3 

2 

<2 

Sensitive  taxa  richness 

>3 

3  -2 

1 

0 

Percent  filterers 

0-5 

5.01  -  10 

10.01- 

25 

>25 

Percent  tolerant  taxa 

0-5 

5.01  -  10 

10.01  - 

35 

>35 

Table  2b.  Provisional  metrics  and  scoring  criteria  for  bioassessment  of  streams  of  the 
Mountain  Ecoregions  (Bukantis  1998). 


Score 

Metric 

3 

2 

1 

0 

Taxa  richness 

>28 

28- 

24 

24  - 

19 

<  19 

EPT  richness 

>  19 

19- 

17 

17- 

15 

<  15 

Biotic  Index 

<3 

3- 

4 

4- 

-5 

>5 

%  dominant  taxon 

<25 

25- 

35 

35- 

45 

>45 

%  collector 

<60 

60- 

70 

70- 

•80       "t 

'             >80 

%  scraper  + 
shredder 

>55 

55- 

40 

40  - 

•  25 

<25 

%EPT 

>70 

70- 

■  55 

55- 

40 

<'4b 

TableSa.  Criteria  for  t±ie  assignment  of  use- support  classifications  /  standards  violation 
thresholds  (Bukantis  1998). 

%  Comparability  to  reference 

Use  support  - 

>75 

25-75 

<25 

FuU  support— standards  not  violated 

Partial  support— moderate  impairment— standards  violated 

Non- support— severe  impaiiiuent— standards  violated 

TableSb.  Criteria  for  the  assignment  of  impairment  classifications  (PlaOdn  et  al.  1989). 

%  Comparability  to  reference 

Classification 

>83 
54-79 
21-50 
<17 

nonimpaired 
slightly  impaired 
moderately  impaired 
severely  impaired 

In  the  narrative  portion  of  this  report,  certain  other  metrics  were  used  as  descriptors 
of  the  benthic  commvinity  response  to  habitat  or  water  quality;  these  metrics  were  not 
incorporated  into  the  revised  bioassessment  index,  either  because  they  have  not  yet 
been  tested  for  reliability  in  stream.s  of  western  Montana,  or  because  results  of  such 
testing  did  not  show  them  to  be  robust  at  distinguishing  impairment,  or  because  they 
did  not  meet  other  requirements  for  inclusion  in  the  mietric  battery.  These  metrics  and 
their  use  in  predicting  the  causes  of  impairment  or  in  describing  its  effects  on  the  biotic 
cormnunity  are  described  below. 

•  The  modified  biotic  index.  This  metric  is  an  adaptation  of  the  Hilsenhoflf  Biotic 
Index  (HBI,  Hilsenhofif  1987),  •\^iich  was  originally  designed  to  indicate  organic 
enrichment  of  waters.  Values  of  this  metric  are  lowest  in  least  impacted 
conditions.  Taxa  tolerant  to  saprobic  conditions  are  also  generally  tolerant  of 
warm  water,  fine  sediment  and  heavy  filamentous  algae  grovrth  (Bollman 
1998b).  Loss  of  canopy  cover  is  often  a  contributor  to  higher  biotic  index  values. 
The  taxa  values  used  in  this  report  are  modified  to  reflect  habitat  and  water 
quality  conditions  in  Montana  (Bukantis  1998).  Ordination  studies  of  the 
benthic  favuia  of  Montana's  foothill  prairie  streams  shovrod  that  there  is  a 
correlation  between  modified  biotic  index  values  and  water  temperature, 
substrate  embeddedness,  and  fine  sediment  (Bollman  1998a).  In  a  study  of 
reference  streams,  the  average  value  of  the  modified  biotic  index  in  least- 
impaired  streams  of  western  Montana  was  2.5  (Wisseman  1992). 

•  Taxa  richness.  This  metric  is  a  simple  count  of  the  number  of  unique  taxa 
present  in  a  sample.  Average  taxa  richness  in  samples  fi-om  reference  streams  in 
western  Montana  was  28  (Wisseman  1992).  Taxa  richness  is  an  expression  of 
biodiversity,  and  generally  decreases  with  degraded  habitat  or  diminished  water 
quality.  However,  taxa  richness  may  show  a  paradoxical  increase  when  mild 
nutrient  enrichment  occurs  in  previousfy  oligotrophic  waters,  so  this  metric 
must  be  interpreted  with  caution. 

•  Percent  predators.  Aquatic  invertebrate  predators  depend  on  a  reliable  source  of 
invertebrate  prey,  and  their  abundance  provides  a  measure  of  the  trophic 
complexity  supported  by  a  site.  Less  disturbed  sites  have  more  plentiful  habitat 

niches  to  suppxjrt  diverse  prey  species,  which  in  turn  suppwrt  abundant  (* 

predator  species. 

o 


•  Nvunber  of  "dinger'  taxa.  So-called  "dinger"  taxa  have  physical  adaptations  that 
allow  them  to  cling  to  smooth  substrates  in  rapidly  flowing  water.  Aquatic 
invertebrate  "dingers"  are  sensitive  to  fine  sediments  that  fill  interstices  between 
substrate  partides  and  eliminate  habitat  complexity.  Animals  that  occupy  the 
hjrporheic  zones  are  included  in  this  group  of  taxa.  Expected  "dinger"  taxa 
richness  in  unimpaired  streams  of  western  Montana  is  at  least  14  (BoUman 
1998b). 

•  Number  of  long-lived  taxa.  Long-lived  or  semivoltine  taxa  require  more  than  a 
year  to  completely  develop,  and  their  numbers  decline  wiien  habitat  and/ or 
water  quality  conditions  are  unstable.  They  may  completely  disappear  if 
channels  are  dewatered  or  if  there  are  periodic  water  temperature  elevations  or 
other  interruptions  to  their  life  cycles.  Western  Montana  streams  with  stable 
habitat  conditions  are  expected  to  support  six  or  more  long-lived  taxa  (Bollman 
1998b). 

RESXJLTS 

Bioassessment 

Figure  2  summarizes  bioassessment  scores  for  aquatic  invertebrate  communities 
sampled  fi-om  sites  on  Bear  Creek,  and  compares  the  results  of  the  2  indices  employed.  ■ 
When  the  Montana  DEQ  index  (Bukantis  1998)  was  applied  to  these  data,  scores 
indicate  that  all  4  studied  sites  on  Bear  Creek  were  slightly  impaired  and  partly 
supported  designated  uses.  Scores  derived  from  the  revised  bioassessment  index 
(Bollman  1998),  on  the  other  hand,  indicated  non-impairment  at  the  3  upstream  sites 
(BEARCOl,  BEAR04,  and  BEAR05)  and  fiall  support  of  uses.  The  lowermost  site, 
BEAR03,  was  moderately  impaired  and  partly  supported  designated  uses.  The  2 
bioassessment  methods  ranked  sites  similarly,  but  scores  differed  by  more  than  10 
percentage  points  for  the  3  upper  sites,  and  by  7  percentage  points  for  site  BEAR03. 

Figure  2.  Comparison  of  total  bioassessment  scores  (reported  as  percent  of  maximuxQ 
score).  The  revised  bioassessment  index  (Bollman  1998a)  was  used  to  calculate  scores 
represented  by  green  bars,  and  the  Montana  DEQ  battery  (Bukantis  1998)  was  used  to 
calculate  scores  represented  by  purple  bars.  Bear  Creek.  2003. 
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Table  4a.    Metric  values,  scores,  and  bioassessments  for  4  sites  on  Bear  Creek,  2003. 
Site  locations  are  given  in  Table  1.  The  revised  bioassessment  method  (Bollioan  1998a) 
was  used  to  calculate  scores. 


SITES 

BEARC 
01 

BEARC 
04 

BEARC 
05 

BEARC 
03 

METRICS 

METRIC  VALUES 

Ephemeroptera  richness 

6 

9 

7 

5 

Plecoptera  richness 

5 

8 

6 

1 

Trichoptera  richness 

8 

5 

5 

10 

Number  of  sensitive  taxa 

3 

4 

4 

1 

%  filterers 

0.94 

0.96 

0.30 

13.51 

%  tolerant  taxa 

1.26 

5.10 

3.64 

15.62 

METRIC  SCORES 

Ephemeroptera  richness 

3 

3 

3 

2 

Plecoptera  richness 

3 

3 

3 

1 

Trichoptera  richness 

3 

3 

3 

3 

Number  of  sensitive  taxa 

2 

3 

3 

1 

%  fUterers 

3 

3 

3 

1 

%  tolerant  taxa 

3 

2 

3 

1 

TOTAL  SCORE  (max.=  18) 

17 

17 

18 

9 

PERCENT  OF  MAX. 

94% 

94% 

100% 

50% 

Impairment  classification* 

NON 

NON 

NON 

MOD 

USE  SUPPORT  t 

FULL 

FULL 

FULL 

PART 

Impairment  classifications:  (NO^f)  non-impaired,  (SLI)  slightly  impaired,  (MOD)  moderately  impaired,  (SEV) 
severely  impaired.  Sec  Table  3b.t  Use  support  designations:  See  Table  3a. 
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Table  4b.   Metric  values,  scores,  and  bioassessments  for  4  sites  on  Bear  Creek.  2003. 
Site  locations  are  given  in  Table  1.  The  Montana  DEQ  bioassessment  method  (Bukantis 
1998)  was  used  to  calculate  scores. 


SITES 

BEARC 
Ol 

BEARC 
04 

BEARC 
05 

BEARC 
03 

METRIC  VALUES 

Taxa  richness 

43 

39 

37 

37 

BPT  richness 

19 

22 

18 

16 

Biotic  index 

3.00 

2.80 

2.86 

4.58 

%  dominant  taxon 

14.78 

33.44 

17.27 

9.91 

%  collector 

46.23 

54.78 

57.88 

55.86 

%EPT 

39.62 

44.27 

52.12 

35.74 

%  scraper  +  shredder 

15.09 

21.02 

23.94 

28.83 

METRIC  SCORES 

Taxa  richness 

3 

3 

3 

3 

EFT  richness 

2 

3 

2 

1 

Biotic  index 

3 

3 

3 

1 

%  dominant  taxon 

3 

2 

3 

3 

%  collector 

3 

3 

3 

3 

%EPT 

0 

1 

1 

0 

%  scraper  +  shredder 

0 

0 

0 

•     1 

TOTAL  SCORE 
(max.=21) 

14 

15 

15 

12 

PERCENT  OF  MAX. 

66.67% 

71.43% 

71.43% 

57.14% 

Impairment 
classification^ 

SLI 

SLI 

SLI 

SLI 

USE  SUPPORT  t 

PART 

PART 

PART 

PART 

ment  classirications:  (NON)  non- 

impaired,  (SLI 

slightly  impai 

red,  (MOD)  mo 

derately  impair 

severely  impaired.  Sec  Table  3b.  t  Use  support  designations:  See  Table  3a. 
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Aquatic  invertebrate  communities 

Interpretations  of  biotic  integrity  in  this  report  are  made  without  reference  to  6- 

results  of  habitat  assessments,  or  any  other  information  about  the  sites  or  watersheds 
that  may  have  accompanied  the  invertebrate  samples.  Interpretations  are  based  entirely 
on:  the  taxonomic  and  functional  composition  of  the  sampled  invertebrate  assemblages; 
the  sensitivities,  tolerances,  physiology,  and  habitus  information  for  individual  taxa 
gleaned  from  the  writer's  research;  the  published  literature,  and  other  expert  sources; 
and  on  the  performance  of  bioassessment  metrics,  described  earlier  in  the  report,  which 
have  been  demonstrated  to  be  useful  tools  for  interpreting  potential  implications  of 
benthic  invertebrate  assemblage  composition. 

The  uppermost  site  on  Bear  Creek  (BEARCO 1)  supported  a  sensitive  and  diverse 
benthic  assemblage.  Six  mayfly  taxa  were  collected  here,  and  the  biotic  index  value 
(3.00)  was  within  expectations  for  a  montane  stream.  These  findings  surest  that  water 
qusdity  was  good  in  this  reach.  Three  cold-stenotherm  taxa  were  present  in  the  sample. 
Cold,  clean  water  appears  to  have  characterized  the  site. 

Eight  caddisfly  taxa  and  18  "dinger"  taxa  were  taken  in  the  sample,  su^esting 
that  sediment  deposition  did  not  appreciably  compromise  stony  substrate  habitats. 
Interstitial  spaces  were  apparently  available  for  colonization;  the  chloroperlid  Paraperla 
sp.,  typically  associated  with  the  h3rporheos,  was  collected.  In  general,  instream 
habitats  were  probably  diverse  and  undisturbed,  since  overall  taxa  richness  (43) 
indicated  high  diversity.  Reach-scade  habitat  features,  such  as  streambanks,  riparian 
zones,  and  channel  morphologic  elements,  were  likely  intact;  the  high  stonefly  taxa 
richness  (5  tax^  may  support  this  hypothesis.  At  least  6  semivoltine  taxa  were  present 
at  the  site,  making  it  seem  unlikely  that  dewatering  was  a  recent  problem  here. 
Scouring  sediment  pulses,  toxic  pollutants,  or  other  catastrophes  that  would  limit  long 
life  cycles  also  appear  unlikely.  Functionally,  the  assemblage  was  skewed  toward 
gatherers,  and  shredders  were  not  as  abundant  as  expected.  Dearth  of  shredders  may 
be  associated  with  limited  riparian  inputs  of  large  organic  material,  or  hydrologic 
conditions  that  do  not  favor  retention  of  such  material. 

Excellent  water  quality  appears  to  have  persisted  farther  downstream,  where  site 
BEARC04  supported  at  least  9  mayfly  taxa  and  the  biotic  index  value  (2.80)  was  low. 
Four  cold-stenotherm  taxa  were  collected  at  this  location,  including  the  stonefly 
Doroneuria  sp.  and  the  mayfly  Drunella  doddsi. 

Sixteen  "dinger'  taxa  and  5  caddisfly  taxa  were  taken  in  the  sample,  indicating 
that  stony  substrates  were  free  of  sediment  deposition.  The  high  taxa  richness  (39)  and 
diverse  predator  fauna  (14  tax^  suggest  that  instream  habitats  of  all  types  were 
abundantly  available.  Eight  stonefly  taxa  were  collected;  high  richness  in  this  group 
may  be  associated  with  stable  streambanks,  intact  riparian  zones,  and  natural  channel 
morphology.  It  seems  unlikely  that  dewatering  or  other  disasters  were  among  recent 
events  in  this  reach,  since  7  long-lived  taxa  were  among  the  sampled  animals.  Here,  as 
at  the  upstream  site,  gatherers  dominated  the  functional  composition  of  the 
assemblage;  shredders  were  only  slightly  more  abundant  at  BEAR04.  The  most 
abundant  (33%)  taxon  at  this  site  was  the  elmid  HeterHmnius  sp.,  a  gatherer.  This 
animal  can  dominate  samples  serendipitously,  since  it  is  gregarious  in  habit. 

At  site  BEAR05,  water  quality  indicators  again  gave  results  suggesting 
unpolluted  conditions  and  appropriate  temperatures.  Seven  mayfly  taxa  were  collected, 
and  the  biotic  index  value  (2.86)  was  characteristic  of  vuiimpaired  montane  streams. 
Three  cold-stenotherm  taxa  were  someu^iat  fewer  than  exp>ected,  but  indicate  that 
water  temperatxires  were  low. 

Five  caddisfly  taxa  and  17  "dinger"  taxa  were  taken  in  the  sample,  suggesting 
that  clean  stony  substrates  free  of  deposited  sediments  were  avaflable.  Thirty-seven 
taxa  in  sdl  appeared  in  the  saiaple,  and  12  of  these  were  predators;  instream  habitats 
were  probably  varied  amd  generally  undisturbed.  Stoneflies  were  weU-represented  (6 
taxa)  which  could  indicate  that  reach-scale  habitat  features  were  essentially  intact.  < 

Stable  streambanks,  functional  riparian  zones,  and  natural  channel  morphological 
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features  may  be  implied.  Six  semivoltine  taxa  were  collected,  including  the  perlid 
stoneflies  Doroneuria  sp.  and  Hesperoperla  pacifica.  Year-round  surface  flow  probably 
characterized  the  reach.  Gatherers  dominated  the  functional  mix,  but  shredders  were 
more  abundant  here  than  at  either  of  the  2  upstream  sites. 

Site  BEARC03  supported  a  benthic  assemblage  that  diEFered  from  those  of  the 
other  Bear  Creek  sites.  Five  mayfly  taxa,  somewhat  fewer  than  expected,  were  collected 
here,  and  the  biotic  index  value  (4.58)  was  elevated  compared  to  typical  values  for 
vmdisturbed  montane  streams.  Tolerant  taxa  made  up  16%  of  the  assemblage.  Not  a 
single  cold-stenotherm  taxon  was  present  in  the  sample,  and  several  taxa  taken  here 
prefer  warmer  water  temperatures.  These  include  the  caddisfly  HeUcopsyche  boreatis 
and  the  mayfly  Tricorythodes  sp.  Water  quality  may  have  been  impaired  by  nutrients  or 
elevated  water  temiperatures.  > 

Ten  caddisfly  taxa  and  1 5  "dinger^  taxa  were  counted,  suggesting  that  stony 
substrates  were  not  obliterated  by  deposited  sediment.  Other  instream  habitats  appear 
to  have  been  undisturbed  since  the  site  supported  at  least  37  distinct  taxa;  nine  of 
these  taxa  were  predators.  A  single  stonefly  taxon,  the  perlodid  Skwala  sp.,  was 
collected.  Low  stonefly  richness  may  be  associated  with  disturbance  to  reach-scale 
habitat  features.  Recent  dewatering  or  other  catastrophes  that  would  abort  long  life 
cycles  seem  ixnlikely,  since  3  semivoltine  taxa  were  taken.  The  functional  composition 
included  all  expected  groups;  the  proportion  of  each  was  probably  appropriate. 

CONCLUSIONS  j 

•  Sites  BEARCOl,  BEARC04,  and  BEARC05  supported  similar  benthic  !  ' 
assemblages.  These  assemblages  were  diverse  and  sensitive,  and  surest  that 
habitat  and  water  quality  conditions  were  good  at  these  sites. 

•  Water  quality  indicators  suggested  that  site  BEARC03  may  have  been  impaired 
by  elevated  water  temperatures  or  mild  nutrient  poUution.  Reach-scale  habitat 
features  may  have  suffered  disturbance  in  this  reach. 

Habitat  Assessment 

Figure  3  graphicaUy  compares  total  habitat  assessment  scores  recorded  for  the  4 
sites  in  this  study.  Table  5  shows  the  habitat  parameters  evaluated,  parameter  scores, 
and  overall  habitat  evaluations  for  the  sites  studied. 

Figure  3.  Total  habitat  assessment  scores  for  sites  on  Bear  Creek,  2003. 
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Table  5.  Stream  and  riparian  habitat  assessment.  Sites  were  assessed  based  upon 
criteria  developed  by  Montana  DEQ  for  streams  with  riffle/ run  prevalence.  Bear  Creek. 
2003. 
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Max. 

possible 
score 

Parameter 

SITES 

BEARC 
01 

BEARC, 
04 

BEARC 
05 

BEARC 
03 

10 

Riffle  development 

8 

6 

9 

9 

10 

Benthic  substrate 

5 

7 

7 

7 

20 

Embeddedness 

5 

14 

16 

15 

20 

Channel  alteration 

15 

20 

18 

20 

20 

Sediment  deposition 

6 

14 

15 

10 

20 

Channel  flow  status 

9 

18 

20 

18 

20 

Bank  stability 

5/5 

6/7 

7/7 

9/10 

20 

Bank  vegetation 

8/6 

8/9 

8/7 

8/10 

20 

Vegetated  zone 

10/6 

5/10 

5/3 

9/10 

160 

Total 

68.00 

101.50 

103.50 

107 

Percent  of  maximum 

52.31% 

78.08% 

79.62% 

82.31% 

CONDITION* 

Marginal 

Optimal 

Optimal 

Optimal 

*  Condition  categories:  Optimal  >  80%  of  maximum  score;  Sub-optimal  75  -  56%;  Marginal  49  -  29%;  Poor 

<23%.  PlafWn  et  al.  1989. 

Overall  habitat  conditions  were  rated  marginal  at  site  BEARCO 1 ,  where 
evaluators  noted  monotonous  benthic  substrate,  embedded  substrate  particles,  and 
moderate  sediment  deposition.  Flow  status  was  rated  marginal  at  this  site. 
Streambanks  vpere  judged  moderately  unstable  on  both  sides  of  the  channel. 

At  BEARC04,  overall  conditions  scored  optimally.  AH  instream  habitat 
parameters  were  perceived  to  be  in  optimal  or  sub-optimal  condition.  The  riparian  zone 
was  abbreviated  or  disturbed  on  the  left  side  of  the  channel. 

Similarly,  instream  habitat  parameters  were  judged  to  be  in  optimal  or  sub- 
optimal  condition  at  site  BEARC05.  The  riparian  zone  width  was  reported  to  be 
marginal  here.  Overall  habitat  conditions  were  rated  optimal. 

At  the  lowermost  site,  BEARC03,  moderate  sediment  deposition  was  reported  by 
evaluators.  Other  instream,  streambank,  and  riparian  zone  habitat  parameters  were 
perceived  to  be  in  optimal  or  sub-optimal  condition  here,  and  the  overall  score  indicated 
optimal  habitat. 

Habitat  assessment  vs.  bioassessment 

When  habitat  assessment  scores  are  plotted  against  bioassessment  scores,  the 
resvdting  figure  provides  an  opportunity  to  evaluate  the  hypothetical  relationship 
between  habitat  integrity  and  water  quality.  Both  factors  are  critical  and  interactive 
determinants  of  the  composition  and  functional  integrity  of  aquatic  invertebrate 
assemblages.  Presumably,  high  quality  habitat,  in  the  absence  of  impairments  to  water 
quality,  supports  functional,  diverse,  and  sensitive  invertebrate  assemblages;  these  are 
assemblages  that  attain  high  bioassessment  scores. 

Barbour  and  Stribling  (1991)  have  hypothesized  that  diminishing  habitat  quality 
should  |>roduce  predictable  diminishment  of  bioassessment  scores,  when  water  quality 
is  not  a  further  insult.  Figure  4  is  a  plot  of  habitat  assessment  scores  against 
bioassessment  scores  (revised  method)  for  the  sampled  assemblages  on  Bear  Creek.  The 
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red  line  superimposed  on  the  plot  represents  the  hypothetical  relationship  between 
habitat  quality  and  biotic  integrity  given  good  water  quality.  In  this  model,  symbols 
falling  in  the  upper  right  area  of  the  graph  would  represent  sites  with  high  scores  for 
both  bioassessment  and  habitat  assessment;  according  to  this  model,  these  would  be 
unimpaired  sites  both  in  terms  of  habitat  integrity  as  well  as  water  quality.  Sites 
BEARCOl,  BEARC04,  and  BEARC05  fall  into  this  region  of  the  plot.  Symbols  falling  on 
or  near  the  red  line  are  hypothesized  to  support  assemblages  affected  by  habitat 
disturbance  in  the  presence  of  adequate  water  quality.  When  habitat  scores  are  high, 
but  bioassessment  scores  are  inordinately  low,  sites  fall  into,  the  lower  right  hand  area 
of  the  plot.  According  to  the  model,  these  sites  support  invertebrate  assemblages  that 
are  impacted  mostly  by  impairment  to  water  quality.  Site  BEARC03  falls  into  this  region 
of  the  graph. 


Figure  4.  Total  bioassessment  scores  plotted  against  habitat  assessment  scores  for 
sites  on  Bear  Creek.  2003.  (Barbour  and  Stribling  1991). 
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